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Abstract 

Tests of local position invariance (LPI) made by comparing the relative redshift of atomic clocks 
based on different atoms have been carried out for a variety of pairs of atomic species. In most 
cases, several absolute frequency measurements per year are used to look for an annual signal, 
resulting in tests that can span on order of a decade. By using the output of continuously running 
clocks, we carry out LPI tests with comparable or higher precision after less than 1.5 years. These 
include new measurements of the difference in redshift anomalies (3 for hyperfine transitions in 
87 Rb and 133 Cs and in 1 H and 133 Cs and a measurement comparing 87 Rb and 1 H, resulting in a 
stringent limit on LPI, /3Rb — /?h = (—2.7 ± 4.9) x 10~ 7 . The method of making these measurements 
for continuous clocks is discussed. 
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I. INTRODUCTION 



Metric theories of gravity are based on the universal coupling of gravity to matter and 
energy, formalized in the Einstein equivalence principle (EEP). This principle states that 
gravitational acceleration is independent of composition (the weak equivalence principle) 
and that nongravitational measurements should be independent of the velocity of the freely- 
falling reference frame (local Lorentz invariance) and of the location in spacetime (local 
position invariance; LPI) where they are carried out. Most efforts to unify general relativity, 
the widely accepted metric theory of gravity, with quantum field theory involve a breakdown 
of EEP |l|,b|. 

According to conventional parametrizations, the principle of LPI is the least-well tested 
tenet of EEP Furthermore, tests of LPI can be used to constrain the coupling of 
certain fundamental constants to gravitational potential and more generally to test for spatial 
dependence of the constants. A spatial dependence is allowed under a particular scenario 



in string theory [2], and there have been published claims of a measured spatial variation of 



the fine structure constant from astronomical observations {4]. 

The traditional redshift for a clock with frequency v depends only on the change in gravi- 
tational potential AU, Avjv = AU/c 2 , and should be independent of the clock composition, 
for example. Violation of LPI may manifest itself in an anomalous gravitational redshift of 
an atomic clock, commonly parameterized with a clock- dependent term 0: 

Avjv = (1 + f3)AU/c 2 . (1) 

Measurements of gravitational redshift vs AU provide a test of the redshift formula and of 
LPI [5|. Alternatively, the relative redshift of two clocks of different composition can be 
used; if the two clocks have LPI- violating parameters (5\ and 02, the relative frequencies 
should vary with gravitational potential as 

(Az//z/) lj2 = (/3i-/32)Af//c 2 . (2) 

This differential measurement can be made with a higher precision than a direct measure- 
ment of redshift vs AU at the possible expense of a reduction in the size of the effect if Pi 
and P2 are close in value. The gain in precision reported here compared to direct redshift 
measurements is about 200, indicating that the two redshift anomalies must be the same to 
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5 x 10 for the best LPI test by absolute measurement of the gravitational redshift to be 



comparable [5] . The possible range of values for fa — fa is one reason why tests with many 
different pairs of clocks is valuable. 

A convenient gravitational potential to use for an LPI test is the solar potential experi- 
enced on earth; due to the earth's elliptical orbit, this varies annually as 

AU s /c 2 = A sm(ut + O ) , (3) 

where A = 1.65 x 10 -10 , u = 0.0172 rad/day, and 0o is the phase such that AU S is a 
minimum at aphelion, such as modified Julian day 55746 (July 4, 2011). (The additional 
variation in potential due to the earth's rotation about its axis is ~ 10~ 3 AU S for our latitude 
and is ignored here.) The difference in LPI- violating parameters for two atomic clocks of 
different composition can be measured by determining the size of the annual oscillation with 
the appropriate phase in the clocks' relative frequency. 

Over the past decade measurements of this sort have been made using a variety of pairs of 
atomic species summarized in Fig. [TJ In almost all previous measurements absolute 

frequencies of a certain type of atomic clock are compared to a cesium standard over many 
years, with typically several (fractional) frequency comparisons at the ~ 10~ 15 level per year. 
We implement a different approach, where each of the atomic clocks used is in continuous 
operation, providing a much higher rate of frequency comparisons and resulting in a more 
efficient LPI test. We measure the output frequencies of clocks based on hyperfine transitions 
in 87 Rb, 133 Cs, and 1 H over 1.5 years. Frequency comparisons between any two of the atomic 
species can be made with a statistical uncertainty of ~ 10~ 15 in a day or two, resulting in 
several hundred data points per year of precision comparable to previous tests. 

II. CLOCKS, METHOD AND ANALYSIS 

Because we do not evaluate absolute frequencies, we have to be confident that our clocks' 
outputs are stable or change in knowable, predictable ways. This is accomplished by main- 
taining a continuous record of a clock's frequency against a variety of other clocks, enabling 
us to determine if discontinuous changes in frequency or drift rate occur. 

The US Naval Observatory (USNO) has a large ensemble of atomic clocks used for precise 
time jl2], including commercial cesium-beam clocks, hydrogen masers, and recently added 
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FIG. 1: (Color online.) LPI tests for various pairs of atomic species in terms of the difference in 
anomalous redshift parameters, A/3 = /3\ — fa- (a) Dark data points are previous measurements: (i) 
neutral strontium optical transition against cesium standard, carried out over 2.5 years [9]fll|: (ii) 
mercury ion optical transition against cesium standard, carried out over 5 years [7]; (iii) hydrogen 
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; and (iv) rubidium-87 



hyperfine transition against cesium standard, carried out over 7 years 

hyperfine transition against cesium standard, carried out over 14 years [u Light data, points are 
results presented in this work, (b) A more detailed look at the shaded region in (a). 



131 ] . All clocks are 



rubidium-fountain clocks, which have been in operation for 1.5 years 
measured regularly against a master reference, which enables a comparison between any two 
clocks. The collection of each type of clock can be used to generate one or more stationary 
outputs representing that ensemble. 

The rubidium frequency record used for LPI tests is an average of the frequencies of 
the two highest performing fountains over the past 1.5 years. Comparisons among rubidium 
fountains demonstrate that these two clocks exhibit extremely stable frequencies with white- 
frequency noise levels below 2 x 10~ 13 . The fountains show no indication of frequency drift 
at the level of 3 x 10~ 18 /day, as determined by using cesium-fountain primary standards 
that contribute to International Atomic Time as a reference. As a baseline, we look for an 
LPI-violating signal in the frequency comparison of the two rubidium fountains. The result 
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is well within the error bar, indicating no bias in measurements using these clocks. 

About 70 commercial cesium clocks are intercompared to create a single output used 
for our LPI analysis [l^. The output of this cesium ensemble is characterized by a white- 
frequency noise level of 1 x 10~ 12 and shows a drift in fractional frequency of 2 x 10~ 17 /day 
over the period in question. There is no long frequency record for an individual cesium 
clock without adjustments based on the performance of the ensemble, so there is no way to 
measure an LPI-violating baseline for the individual clocks. 

Hydrogen masers generally exhibit complicated behavior over long averaging times 
(months). They are known to be sensitive to environmental factors, which can cause fre- 
quency discontinuities in addition to less violent perturbations, and they exhibit significant 
frequency drifts, which may not be linear and which can change in time. Of the ~ 25 
hydrogen masers available to us, a subset is ruled out for these LPI tests because they ex- 
hibit discontinuous changes in frequency or drift rate during the interval of time that we 
are interested in. Further assessing which masers to include requires careful analysis of the 
frequency records of different masers compared against each other as well as against the 
rubidium and cesium clocks used in the analysis, as discussed further below. The maser 
short-term frequency stability is better than that of either of the other types of clock used 
for these LPI tests. The white-frequency noise level of a maser-fountain comparison is typ- 
ically of the order of 5 x 10~ 13 , limited by noise in the frequency-counter measurement 
system. Some masers are also measured with a dual-mixer measurement system, in which 
case a maser-fountain comparison exhibits a white-frequency noise level limited by fountain 
performance. 

In order to measure /?i — 02 for a pair of atomic species we need to examine the relative 
frequency record for an annual oscillation with phase <po in the presence of a drift. We 
limit our analysis to clocks that show a frequency drift that is linear, which, for a ~ 1.5- 
year interval, has little impact on the determination of the sinusoid amplitude. Applying a 
least-squares fit using a function that is a combination of drift and oscillation gives results 
identical to fitting the frequency record to a line and then fitting those residuals to a sinusoid. 
We use the latter procedure for presenting data in the figures below. All relative frequency 
measurements participating in the fit derive from averaging intervals resulting in white 
frequency noise, consistent with application of least-squares curve fitting. 
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FIG. 2: (Color online.) (a) Frequency difference of rubidium-fountain average and cesium ensemble 
over 1.5 years. Fit to a linear drift is shown, (b) Residuals after removing drift, along with the fit 
to an LPI-violating signal, amplified x20. MJD stands for modified Julian day. 

A. Rubidium vs Cesium 

Application of this procedure to a rubidium-cesium LPI test is straightforward. Fig- 
ure EJa) shows the fractional frequency difference between the cesium ensemble and the 
average of rubidium fountains over 1.5 years. Raw data are recorded at hourly intervals; 
the plots in the figure show the frequency averaged over 2 days. The graph includes a fit 
indicating a relative drift of 2 x 10~ 17 /day. In Fig. [2]^b), the residuals after removing the 
drift are fit to an LPI-violating signal, shown amplified in the graph. The amplitude for this 
fit is Ai/Rb/^Rb — A^Cs/Vcs = (2.6 ± 2.1) x 1CT 16 (all uncertainties are 1 standard deviation). 
This translates to a difference in redshift anomalies of 

/3 R b-/3cs = (-1.6±1.3)xl(r 6 . (4) 

B. Hydrogen 

From the subset of masers that exhibit no apparent discontinuities in frequency or drift 
over all or a large portion of the past 1.5 years, we need to establish which masers have 
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frequency records that are truly linear, without higher order terms. Some masers obviously 
exhibit frequencies that change in time in a nonlinear way and can be removed from consid- 
eration. To rule out masers with a more subtle nonlinear behavior, we rely on the quality 
of the least-squares fit to the relative frequency record of a maser and the rubidium aver- 
age [15]. If the data do not fit well, as indicated by the reduced x 2 an d a visual inspection of 
the residuals, the maser is not included in the analysis. We attribute a poor fit to unmodeled 
behavior in the maser; this procedure leaves us with four masers in the LPI analysis. 

It is possible that these four masers still display unmodeled behavior, which mimics an 
LPI-violating signal. Such unmodeled behavior would not be surprising since the phase of 
the LPI-violating signal is close to that of the seasonal cycle. We asses this by looking at all 
six possible pairs of the four masers for an LPI-violating signal. These measurements enable 
a (3 to be determined for each maser that serves as a baseline for remaining unmodeled 
behavior. Each value serves as a bias and is subtracted from the measured LPI-violating 
signal for each maser-rubidium or maser-cesium comparison. Once this bias is included in 
a clock comparison using a given maser, the final uncertainty for that comparison is the 
combined statistical uncertainty from the LPI fit and the uncertainty associated with the 
bias. 

In order to account for the different uncertainties associated with each maser, the hydro- 
gen results are determined by measuring LPI-violating amplitudes for each of the four masers 



against the cesium and rubidium clocks and combining them in a weighted average 161 ] . The 
four measurements against the cesium ensemble give 

/3h-/3 Cs = (-0.7±1.1)x10- 6 . (5) 

The average of the rubidium fountains measured against the four masers yields 

/3Rb-/3 H = (-2.7±4.9)xl0- 7 . (6) 

Figure [3] shows plots of the frequency of one of the masers compared against [Fig. [31(a)] the 
cesium ensemble and [Fig. [2(b)] the rubidium average. 



III. DISCUSSION AND SUMMARY 



Measurements of fii — $2 for different pairs of atomic species can be used to constrain 
the coupling to gravitational potential of certain dimensionless constants: the fine structure 
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FIG. 3: (Color online.) Frequency comparison for one hydrogen maser and (a) the cesium ensemble 
and (b) the rubidium-fountain average. Each point is a 2-day average, and linear drifts have been 
removed. Fits to LPI-violating signals are amplified by the factors shown. 

constant a, the electron-to-proton mass ratio m e /m p , and the ratio of the light quark mass 
to the quantum chromodynamics length scale, m q /AqcB, where m q is the average of the up 
and down quark masses 17| . For the case of massive bodies the size of the sun, the coupling 
of these constants to gravity is expected to be proportional to the coupling to a hypothetical 
scalar field that is a component of many cosmological models and that could affect their 
values 

The coupling of these fundamental constants to gravitational potential can be character- 
ized by dimensionless coupling parameters k e fl?) ]: 

/ACA 



e 



k e 



V c 2 



(7) 



where e stands for the three constants a, m e /m p and m q /AQcr>- Using calculated sensitiv- 
ities of different atomic transition frequencies to variations in these constants, differences 



in redshift anomalies for a pair o 
the coupling parameters k. 
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atomic species can be related to linear combinations of 



19| . The three LPI measurements presented here can be 



used to put limits on the coupling of a and m^/AgcD using only clocks from one institu- 
tion, k a = (3.3 ±2.3) x 10~ 6 and k mq /^ QCD = (—1.8 ± 1.3) x 10 -5 . Measurements between 
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FIG. 4: (Color online.) Constraints on coupling of a, m e /m p and m g /AQCD to gravitational 
potential characterized by the dimensionless parameters k a , k me i m , and L a . Gray bands 
show constraints imposed by previous measurements, compiled from Refs. 7H10l| and presented 
in [l(J. Diamonds show the tighter constraints that result when our measurements are included 
in the analysis. Squares represent the values for k a and &m 9 /A QCD that are derived using only the 
results presented here. 

optical and microwave clocks are needed to constrain the coupling for m e /m p . By using our 
measurements along with those from references 7H10|| we get limits on all three couplings: 



k 



(1.7 ±7.5) x 1(T 7 
(-2.5 ± 5.4) x 1(T 6 
(3.8 ±4.9) x 1(T 6 . 



(8) 



These results are plotted in Fig. HJ 

The tests involving cesium are limited by the white-frequency noise of the cesium ensem- 
ble, so the precision of these measurements should improve with longer fitting times. The 
limit on the precision of the rubidium-hydrogen test is a combination of white-frequency 
noise from the measurement system and uncertainty in the LPI-mimicking amplitude for 
each maser. While these uncertainties should also improve with longer fitting times, the 
unpredictable nature of masers makes this less certain. Even now, with an uncertainty on 
the rubidium- hydrogen measurement in the mid 10~ 7 s, this very precise LPI test could re- 
main competitive for the near-future. While relative frequency measurements of two optical 
clocks should yield much more stringent limits, many optical clock frequencies are mea- 
sured against the cesium standard, limiting the precision to a level governed by a microwave 



transition. Furthermore, our continuous clock measurements are competitive with LPI tests 
that could be implemented in near-term space missions. Clock experiments intended for 



the International Space Station [20J could improve upon the precision of absolute redshift 
measurements but not differential measurements; the driving term due to the solar potential 
is the same as on earth, and there is no driving term from the earth's potential due to the 
station's circular orbit. In terms of space clocks, improvements to differential redshift tests 



of LP 
orbit 



could result from proposed missions that would put clocks in a highly eccentric earth 



2l| or a solar system escape trajectory [22]. 



In summary, we have made precise LPI tests with three different pairs of atomic species 
using continuously running atomic clocks and have used these to tighten the constraints on 
coupling of fundamental constants to gravity. 
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